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HEAT  AND  MASS  TRANSFER  FROM 
FREELY  FALLING  DROPS  AT  LOW 
TEMPERATURES 


John  P.  Zarling 

INTRODUCTION 

Icc  has  been  used  numerous  times  as  a  structural  material  in  cold  regions  to  substantially  reduce 
construction  costs  for  roads,  bridges,  and  work  pads.  To  shorten  construction  time,  techniques 
such  as  surface  flooding  or  spraying  of  water  have  been  used  in  the  past,  as  reported  by  Hoffman 
(1967)  and  McFadden  (1976),  to  accelerate  ice  growth.  However,  very  little  engineering  design  in¬ 
formation  exists  related  to  ice  growth  rates  that  can  be  achieved  using  these  artificial  thickening 
methods.  The  rapid  ice  building  project  currently  being  conducted  by  the  Alaskan  Projects  Office 
of  CRREL  is  an  attempt  to  provide  an  engineering  design  methodology  which  will  predict  the  rates 
of  artificial  icc  thickening. 

The  approach  being  taken  by  the  Alaskan  Projects  Office  is  two-fold:  experimental  and  analyti¬ 
cal.  A  pump-nozzle  system  (sec  cover  photo)  has  been  designed,  assembled,  and  used  to  pump  0°C 
water  from  beneath  the  icc  sheet  and  jet  it  into  the  air.  As  the  water  jet  breaks  up  into  drops,  a 
large  surface  arca/volume  ratio  results  which  produces  high  heat  transfer  rates  and  subsequent  super¬ 
cooling  of  the  drops  as  they  fall  back  to  the  ice  surface.  In  fact,  field  tests  have  shown  that  some 
drops  actually  nucleate  into  icc  before  striking  the  ice  sheet.  This  report,  however,  will  only  exam¬ 
ine  the  heat  and  mass  transfer  rates  from  freely  falling  drops  on  an  analytical  engineering  design  basis. 

There  have  been  several  research  programs,  both  theoretical  and  experimental,  conducted  on  heal 
and  mass  transfer  from  drops.  The  bulk  of  these  studies  have  cooling  tower,  spray  cooling,  or  chem¬ 
ical  processing  applications  as  their  focus.  Ranz  and  Marshall  (1952)  experimentally  investigated 
the  evaporation  rates  from  drops,  their  studies  confirming  the  analogy  between  heat  and  mass  trans¬ 
fer.  They  also  presented  correlative  expressions  for  the  Nusselt  and  Sherwood  numbers  in  terms  of 
Reynolds  and  Prandtl  or  Schmidt  numbers.  Yao  and  Schrock  (1976)  conducted  a  scries  of  experi¬ 
ments  on  heat  and  mass  transfer  from  freely  falling  drops  over  the  zero  to  terminal  velocity  range. 
They  modified  the  Ranz  and  Marshall  expressions  to  include  acceleration  effects.  Chao  (1969)  per¬ 
formed  an  analytical  study  of  transient  heat  and  mass  transfer  from  translating  drops  at  large  Rey¬ 
nolds  and  Peclet  numbers.  He  included  internal  circulation  within  the  drop  and  performed  a  bound¬ 
ary  layer  similarity  analysis  at  the  drop  surface  to  arrive  at  transient  and  steady-state  Nusselt  numbers. 

The  behavior  of  drops  and  sprays  from  a  d\  namic  point  of  view  has  also  received  attention.  Hughes 
and  Gilliland  (1952)  presented  a  general  review  including  the  gross  motion  of  drops  and  the  detailed 
motion  in  and  around  individual  drops.  The  emphasis  of  the  article  is  on  correlations  in  connection 
with  the  effect  of  acceleration  on  drag,  the  equilibrium  distortion  or  oscillations  from  prolate  to 
oblate  spheroidal  shapes,  and  the  internal  circulation  caused  by  skin  friction.  Shattering  of  drops 
due  to  aerodynamic  forces  under  both  transient  and  steady  state  conditions  was  studied  experiment¬ 
ally  by  Lane  (1951).  The  largest  drops  he  tested  were  4  mm  in  diameter  which  required  a  steady 
relative  air  speed  of  1 2  m/s  before  shattering  occurred.  Smaller  drops  required  higher  shattering 
velocities.  (The  free  fall  speed  of  4-mm-diam  water  drops  is  less  than  the  shattering  velocity.)  Lapplc 
and  Sheperd  (1940)  have  presented  equations  and  drag  coefficients  for  calculating  the  trajectories  of 
individual  spherical  particles  undergoing  accelerated  motion  due  to  gravity  in  a  viscous  fluid.  Exper¬ 
imental  tests  also  conducted  by  them  compared  the  actual  and  theoretical  trajectories,  which  were 
in  fairly  good  agreement. 


Figure  1.  Nucieation  temperature  of 
water  drops  as  a  function  of  water 
purity  and  drop  diameter. 


The  rate  of  evaporation  from  sprays  has  been  reported  by  Dickinson  and  Marshall  (1 968).  The 
principal  parameters  considered  were  the  drop  mean  diameter  and  size  distribution,  initial  drop  vel¬ 
ocity  and  temperature,  and  air  velocity  and  temperature.  Computations  were  carried  out  and  results 
reported  on  the  interrelationship  of  these  parameters.  Chen  and  Trczck  (1977)  introduced  a  wet 
bulb  temperature  weighting  factor  which  was  applied  to  the  single  drop  thermal  performance  model 
to  account  for  local  temperature  and  humidity  variations  within  the  spray  field. 

Nucieation  and  dendritic  ice  formation  within  water  drops  falling  through  low  temperature  air  are 
related  to  the  formation  of  icc  fog  produced  by  water  vapor  output  from  automobile  exhaust,  power 
plant  stacks,  residential  chimneys,  power  plant  cooling  ponds  and  other  sources  associated  with  an 
urban  environment.  Vapor  fog  from  these  sources  nucleates  to  an  icc  fog  whenever  the  ambient 
temperature  approaches  -30° C.  Benson  (1970)  presented  an  overall  review  of  the  icc  fog  phenome¬ 
non  including  its  sources,  distribution,  and  physical  nature  in  interior  Alaska.  A  section  of  his  report 
discusses  the  supercooling  and  nucieation  of  small  water  drops  into  ice.  An  extensive  review  of  nu- 
clcation  temperatures  of  supercooled  water  drops  has  been  conducted  by  Langham  and  Mason  (1958). 
Figure  1  shows  that  the  nucieation  temperature  of  water  drops  decreases  with  decreasing  drop  size 
and  increasing  water  purity  until  the  homogeneous  nucieation  temperature,  -40°C,  is  reached.  Once 
nucieation  has  occurred,  the  drop  temperature  rises  to  0°C  as  the  latent  heat  of  fusion  liberated  dur¬ 
ing  ice  formation  warms  the  supercooled  water.  This  results  in  an  increased  heat  transfer  rate  after 
nucieation  due  to  the  increased  temperature  difference. 

Another  application  of  this  analysis  is  the  purification  of  salt  water  (sea  water)  by  spraying  water 
into  the  air  under  winter  arctic  conditions  (Peyton  and  Johnson  1967).  The  ice  dendrites  formed 
following  nucieation  of  the  falling  drops  are  fresher  or  less  saline  than  the  salt  water.  After  the  ice- 
water  drops  land,  the  brine  drains  away,  leaving  a  mass  of  low  salinity  ice  particles  that  might  be 
used  as  a  domestic  water  source. 

The  remainder  of  this  report  has  been  sub-divided  into  three  sections:  the  velocity  problem,  heat 
and  mass  transfer  from  a  single  drop,  and  heat  and  mass  from  a  system  of  drops.  The  solution  to  the 
velocity  problem  is  required  for  the  calculation  of  the  heat  and  mass  transfer  rates  for  both  the  single 
drop  and  system  of  drops.  The  single  drop  solution  provides  insight  into  the  maximum  rates  of 
heat  transfer  and  mass  transfer.  Flowever,  this  model  does  not  account  for  the  changes  in  the  air 
temperature  and  humidity  as  the  drops  fall.  These  effects  have  been  included  in  the  model  of  a  sys¬ 
tem  of  drops  which  results  in  a  decrease  in  the  cooling  rate  of  the  drops  as  the  ambient  air  tempera¬ 
ture  and  humidity  surrounding  the  drops  increase.  In  practice,  spraying  water  with  a  sweeping  mo¬ 
tion  of  the  nozzle  will  aid  in  enhancing  the  cooling  rate  from  the  drops  as  compared  to  a  fixed  nozzle 
arrangement,  especially  during  calm  wind  conditions. 
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VELOCITY  PROBLEM 


In  order  to  predict  the  heat  and  mass  transfer  rates  from  the  water  drops,  the  velocity  of  the  drops 
must  be  determined.  The  water  jet  trajectory  obeys  the  laws  of  Newtonian  dynamics.  However,  the 
problem  is  complicated  by  the  fact  that  the  jet  breaks  up  into  individual  drops  near  the  apex  of  its 
path  due  to  the  interactions  of  inertial,  aerodynamic  and  surface  tension  forces.  In  the  horizontal 
direction,  the  speed  of  the  drops  will  decrease  continuously  from  the  point  of  discharge  at  the  nozzle 
to  the  point  where  the  drops  strike  the  ice  due  to  aerodynamic  drag.  In  the  vertical  direction,  the  jet 
will  have  its  maximum  speed  at  the  nozzle,  decrease  to  zero  velocity  at  the  apex  of  the  trajectory 
where  the  jet  breaks  up,  and  then  the  drops  will  increase  in  speed  to  terminal  velocity  as  they  fall 
downward. 

The  equation  of  motion  for  an  individual  drop  of  diameter  d  in  vector  notation  (Lapple  and  Shep- 
erd  1940)  can  be  written  as 


dU  +  3 
dr  4  d 


=  0 


0) 


where  U  =  velocity 
r  =  time 


Cp  =  drag  coefficient  that  is  a  function  of  its  Reynolds  number 
d  =  drop  diameter 
g  =  gravitational  acceleration 
pa  =  air  density 
pw  =  water  drop  density. 


At  terminal  velocity  in  they  direction,  dUJdr  is  zero  and  the  resulting  equation  can  be  solved  for 
Uy  as 


Uy  = 


dg(pw-pa) 


3PaCD 


(2) 


For  liquid  droplets,  the  following  drag  coefficient  correlations  are  available  (Chen  and  Trezek  1977), 


24  +  6 

Re 


0.27  1  <  Re  <  1000 


(3) 


and 


CD  =  0.6649 -0.271  2  x  10"3Re  +  1 .22  x  10'7Re2-10.919x  10"nRe3 


(4) 


for  1000  <  Re  <  3600,  where  Re  is  the  Reynolds  number. 


Based  on  eq  2  and  3  or  4,  terminal  velocities  of 
drops  of  varying  diameter  falling  vertically  have  been 
calculated  at  an  air  temperature  of  -1 8°C  and  are  listed 
in  Table  1. 

In  all  subsequent  calculations  in  this  report,  the  vel¬ 
ocity  of  the  freely  falling  drops  is  assumed  to  be  equal 
to  their  terminal  speed  in  the  vertical  ory  direction. 
Reynolds  numbers  will  also  be  based  on  the  terminal 
velocity.  Field  observations  of  jetting  water  into  the 
air  supports  this  assumption,  since  the  horizontal  vel¬ 
ocity  component  of  the  drops  appears  very  small  be¬ 
yond  the  apex  of  the  jet,  especially  at  large  vertical 
angles. 


Table  1.  Terminal  speed  vs  drop  diameter. 

D  U 


Drop 

diameter 

(mm) 

cD 

Drag 

coefficient 

Terminal 

velocity 

(m/s) 

Re 

Reynolds 

number 

0.5 

1.129 

2.07 

87 

1.27 

0.583 

4.54 

477 

2.54 

0.495 

6.95 

1463 

3.81 

0.577 

7,86 

2489 

5.08 

0.727 

8.1 1 

3417 

6.35 

0.889 

8.20 

4317 
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Figure  2.  Heat  and  mass  transfer  be¬ 
tween  water  drop  and  ambient  air. 


HEAT  AND  MASS  TRANSFER  PROBLEM-A  SINGLE  DROP 

There  are  two  approaches  which  have  been  used  in  the  past  to  analyze  the  heat  and  mass  transfer 
rates  from  a  freely  falling  drop.  The  first  approach  is  to  assume  that  the  internal  motion  within  the 
drop  is  so  vigorous  that  complete  mixing  occurs.  The  temperature  profile  is  flat  and  the  thermal 
and  mass  transfer  resistances  occur  only  at  the  drop  surface.  The  alternate  approach  is  to  assume  no 
mixing  in  the  drop  so  that  the  energy  equation  reduces  to  a  transient  heat  flow  problem.  Intuitively, 
the  mixing  model  will  yield  answers  indicating  a  higher  heat  flux  than  the  nonmixing  model. 

The  governing  differential  equation  for  the  heat  and  mass  transfer  from  a  falling,  well-mixed  drop, 
as  shown  in  Figure  2,  is 

PwCp  V  d-~  =  -A\hc( fw-  7"J  +/tdpa/?fg(lfr  WJ  +  oe( T*-T*)\  (5) 


where  Cp 
V 
A 

hc 

K 

h<i 

Pa 

Pw 


nH 

M' 


a 

e 


specific  heat  of  water 

drop  volume 

surface  area  of  drop 

a  convective  heat  transfer  coefficient 

water  drop  temperature 

ambient  temperature 

a  convective  mass  transfer  coefficient 

air  density 

water  density 

enthalpy  of  vaporization  of  water 
saturation  humidity  ratio 
ambient  humidity  ratio 
Stefan-Boltzmann  constant 
emissivity. 


The  term  on  the  left-hand  side  of  the  above  equation  represents  the  time  rate  of  change  of  energy 
of  the  drop  as  it  falls.  The  three  terms  on  the  right-hand  side  of  this  equation  represent  thermal 
energy  transports  by  convection,  evaporation  and  radiation,  respectively. 

In  order  to  solve  this  differential  equation,  the  heat  and  mass  transfer  coefficients  will  be  assumed 
constant.  Furthermore,  the  fluid  physical  properties  will  also  be  assumed  invariant.  The  radiation 
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.cun  in  the  governing  equation  can  ho  simplified  In  introducing  the  radiation  heat  tiansloi  coellici- 
enl  which  is  dolinod  as 


h,  ^  oe(  /,;  +  £)(/„  <  /  ,  ).  (6) 

I  or  small  temperature  diflerences,  /?,  can  he  approximated  as  a  constant  or 

lir  =  oe(r£  +  lJ)(  Ta  +  /„ )  (7) 

where  Ta  is  the  initial  drop  temperature.  An  analysis  ol  Mcl  addon's  (1976)  data  on  arctic  cooling 
pond  water  surface  temperatures,  air  temperatures,  and  net  outgoing  radiation  indicates  that  the 
ambient  aii  temperature  closely  approximates  the  effective  sky  tempciature.  Therefore,  the  ambient 
air  tempciature  was  chosen  as  the  ellective  skv  temperature  to  which  the  drop  radiates. 

Next  the  evaporation  term  must  he  expressed  in  terms  ol  temperature.  Since  the  saturation  humid¬ 
ity  ratio  is  a  function  ol  temperature,  a  parabolic  curve  using  temperature  as  the  independent  variable- 
can  he  used  to  lit  humidity  ratio  data  at  saturation  ovei  the  appropriate  temperature  range, 


Ws-a7$+bT„  +  c,  (8) 

or  in  terms  of  the  humidity  ratio  potential, 

K’w-  ~a(Tj-72)+b(T„-rj.  (9) 

The  coefficients  in  the  above  equations  have  been  determined  over  three  temperature  ranges  using  a 
least-squares  fit  of  the  humidity  ratio  at  saturation  conditions  and  are  given  below: 

ft;  -  6 .094 8 x  I0'ft  T-  +  2.7723 x  I  O'4  rw  +  3.7653x  10'  ’  -20UC  „  T  0”C  +  4ua 

K's=  5.I896x  lO'6  r*  +  2.61 84x  I0'4TVV  +  3.7382x10'-'  -25°C  *.  T  «,  0°C  i  8% 

fks  =  4.4299 x  10'6  7;  +  2.4589x  I0'4  7W  +  3.7002x  10"-’  -30°C  ^  7  s  0°C  ±14%. 


It  a  more  accurate  representation  ol  the  humidity  ratio-temper alure  (unction  is  desired,  a  higher 
order  polynomial  can  be  used  or  the  quadradic  equation  can  be  "patched" over  smaller  temperature 
intervals. 

There  are  two  approaches  which  can  be  used  to  arrive  at  the  mass  transfer  coefficient.  The  lirst 
uses  the  Lewis  number,  Le,  to  relate  the  forced  convective  heat  and  mass  transler  coefficients  for 
water  in  air  (Threlkeld  1970),  or 


Le  = 


1  a'-'p.a^d 


(10) 


where  a  is  the  thermal  ditfusivity  and  D  the  mass  diffusivjty.  For  a  water-air  system,  a/D  is  approxi¬ 
mately  0.85,  from  which  the  value  of  the  convective  mass  transfer  coefficient  hd  can  be  easily  cal¬ 
culated  from  known  correlations  for  the  convective  heat  transfer  coefficient.  The  second  and  seem¬ 
ingly  more  direct  approach  is  to  use  the  Ranz  and  Marshall  (1952)  convective  heat  and  mass  transfer 
coefficients  for  drops  which  have  been  correlated  from  experimental  data  as 


Nu  =  ^  =  2  +  0.6Prl/3  Re"2 
k 


(11) 


and 

5h  =—  =  2+ 0.65c1  3  Re1  1 
D 

5 


r  * 


(12) 


where  Nu  -  Nusscll  number 
Pi  -  Prandtl  number 
Sh  =  Sherwood  number 
Sc  =  Schmidt  number. 


Lquation  5  can  now  be  rewritten  as 


i7w-^)+/,dPa/'fg°(7w2-ri 

OT  Pw'-n« 


(13) 


w  p 


To  simplify  the  algebra  in  the  above  equation,  let 


p  =  6<*c**r *Wn*)  .  c=6pA1v  0  =  T  -T 


PwCprf 


PwCp^ 


After  introducing  these  definitions  into  cq  13  and  rearranging,  the  following  equation  results; 


^  =  -B0-C0(0+2TJ. 
dr 


(14) 


Integrating  the  above  equation  and  then  applying  the  initial  condition  that  0  =  0 j(  (Tw  =  7)),  at 
r  =  0  yields  the  solution 


0  = 


70, 


yeyr+C0,(Uc^) 


(15) 


where  7  =  fi+27"„C. 

The  amount  of  supercooling  of  drops  initially  at  0°C  ranging  in  size  from  0.5  mm  to  6.5  mm  and 
falling  through  ambient  air  at  -1 8°C  was  calculated  based  on  cq  1 5  with  the  results  shown  in  Figure 
3.  As  expected,  the  smaller  drops  cool  to  a  lower  temperature  when  falling  a  given  distance  as  com¬ 
pared  to  larger  drops.  This  fact  is  chiefly  due  to  the  larger  volume/surface  area  ratio  of  larger  drops 
and  the  increased  falling  time  of  small  drops.  The  ratio  of  energy  lost  by  the  drop  to  the  latent  heat 
of  fusion  is  easily  estimated  as  Cp(7'f-  T0)/L  where  Tf,  TD and  L  are  the  freezing  temperature,  final 
drop  temperature,  and  latent  heat  of  fusion,  respectively. 

The  relative  amounts  of  cooling  by  convection,  radiation  and  evaporation  from  the  drop  can  be 
determined  by  evaluation  of  these  respective  terms  in  eq  5.  For  a  5-mm-diam  drop  falling  at  ter¬ 
minal  speed  at  0°C  in  a  saturated  -1 8°C  ambient  air  environment,  convective  heat  transfer  is  68%, 
radiative  heat  transfer  is  2%  and  evaporative  heat  transfer  is  30%. 

In  order  to  compare  the  effectiveness  of  jetting  water  into  the  air  to  surface  flooding,  the  heat 
transfer  rates  of  the  two  techniques  must  be  evaluated.  The  heat  loss  for  the  surface  flooding  case 


Figure  3.  Cooling  curves  for  water 
crops  falling  from  various  heights  in 
a  -I8°C  environment. 
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can  be  determined  by  following  MeFadden's  (1976)  recommendations  on  heat  transfer  equations 
for  calculating  heat  loss  from  cooling  ponds  in  low  temperature  environments.  The  Rimsha-Don- 
chcnko  equation  for  predicting  the  free  convective  heat  loss  from  open  water  during  winter  is  given  as 

Qc  =  [3.9+o.i7(rw-ra)|(rw-ra)  (w/m2).  (16) 

Dalton’s  equation  in  modified  form  to  account  for  evaporative  heat  loss  at  low  ambient  temperatures 
is  presented  as 

Qe  =  21 ,820 (W^-  Wa)  (W/m2).  (17) 

The  radiative  loss  can  be  predicted  in  an  identical  manner  as  was  done  for  the  drops,  or 

Q,  =  oe(T4-T4)  (W/m2).  (18) 

Assuming  no  wind,  a  water  temperature  of  0°C,  and  an  ambient  air  temperature  of  18°C,  the 
water  surface  heat  losses  by  convection,  evaporation,  and  radiation  are  1 25  W/m2,  65  W/m2,  and 
89  W/m2,  respectively.  For  a  4-mm-diam  drop  at  0°C  falling  at  terminal  speed,  convective,  evapora¬ 
tive,  and  radiative  heat  loss  rates  are  1024  W/m2, 375  W/m2,  and  89  W/m2,  respectively.  Added  to 
this  is  the  fact  that  the  drops  provide  a  much  larger  exposed  surface  area  per  unit  volume  of  water 
than  does  ponded  water,  which  greatly  enhances  the  cooling  of  the  sprayed  water  and  subsequently 
increases  the  ice  formation  rate.  For  example,  a  5-cm  layer  of  water  on  the  ice  surface  would  con¬ 
tain  the  equivalent  of  approximately  1 .5  million  4-mm-diam  drops  per  meter2  of  surface  area.  This 
yields  a  drop  surface  to  water  surface  area  ratio  of  75  to  1 . 

HEAT  AND  MASS  TRANSFER-A  SYSTEM  OF  DROPS 

Figure  4  shows  a  schematic  diagram  of  the  column  of  cold  air  and  a  falling  water  spray.  Assum¬ 
ing  steady-flow  conditions,  an  energy  balance  on  a  differential  volume  element  yields 

madh  =  -  +dW-  W-^\  dhf+mahfdW  (19) 

which  can  be  simplified  to 

madh  = -mvdhf+mahfdW  (20) 


Figure  4.  Schematic  diagram  of  energy-mass  balance  across  the  air  column  and  differ¬ 
ential  element. 
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by  neglecting  the  second  order  terms  in  eq  1 9.  In  eq  20,  rni  is  the  mass  flow  rate  of  air,  mw  is  the 
mass  flow  rate  of  water,  dh  is  the  change  in  enthalpy  of  the  air,  dh f  the  change  in  enthalpy  of  the 
water,  and  dW  the  change  in  humidity  ratio  of  the  air  across  the  differential  volume.  The  change  in 
energy  of  the  water  in  the  control  volume  is  due  to  convective  heat  and  mass  transfer  from  the 
drops  or 


-m^dhf=hcAwdV(Tw-Ti)  +  pahdAvdV(Wi-W)hf6.  (21) 

The  term  on  the  left-hand  side  of  the  above  equation  represents  the  change  in  energy  of  the  water, 
and  the  two  terms  on  the  right-hand  side  of  the  equation  represent  the  heat  transfer  due  to  convec¬ 
tion  and  evaporation,  respectively.  In  this  equation,  hc  and  hd  are  the  heat  and  mass  transfer  coeffi¬ 
cients,  Av  is  the  surface  area  of  water  drops  contained  in  the  control  volume,  Tw  and  Ta  are  the  water 
drop  and  air  dry  bulb  temperatures,  and  / 7fg  is  the  enthalpy  of  evaporation  of  the  water  spray.  The 
change  in  concentration  of  water  vapor  in  air  must  be  equal  to  the  mass  transfer  from  the  drops  or 

m1dW  =  pahdAvdV(Wi-W)  (22) 

where  If,  and  W  are  the  humidity  ratios  of  the  saturated  and  moist  air  streams.  Upon  substituting 
the  Lewis  number,  Le  =  /tc/hd  paCp  a  into  eq  21 ,  the  following  equation  is  obtained 


/ 


-mwd/7f  pahd A vd V  [  LeCp^a 


(Tw-Ta)  +  (WS-W)hfg\. 


Substituting  eq  20  and  22  into  eq  23  and  rearranging  yields 


dh_ 

dW 


LeCP,a 


iK-Ta) 

(Ws-W) 


+v 


(23) 


(24) 


Now  using  the  approximation  of  constant  specific  heat  Cp  a  in  the  defining  equation  for  the  enthalpy 
of  the  air  stream,  or 


h  =  Cp  i  Ta  +  2501  W 


eq  25  can  then  be  used  to  write 


(25) 


CpJTw-Ta)  =  (h,-h)  -  2501  {W^W).  (26) 

Equation  26  is  now  substituted  into  eq  24  to  yield 

d±  =  Le  1-lH  *  (/?_  -  2501  Le) .  (27) 

dW  (Ws-W)  g 

This  equation  represents  the  operating  or  condition  line  on  the  psychrometric  chart  for  the  changes 
in  state  as  water  drops  pass  through  the  column  of  air.  Although  data  for  Le  for  this  application  are 
not  available,  Threlkeld  (1970)  recommends  an  approximate  relationship  for  air  flow  given  previously 
by  eq  10.  Thermal  and  vapor  diffusivity  data  for  dry  and  saturated  moist  air  can  be  found  in  Eckert 
(1959).  An  accurate  although  approximate  solution  to  eq  27  can  be  obtained  graphically  on  a  psy¬ 
chrometric  chart.  For  one  set  of  operating  conditions  for  the  air  column,  assume  that  the  entering 
water  temperatures  rwl,  the  water  flow  rate  mw,  the  air  flow  rate  mv  and  the  initial  air  temperature 
Tao  and  pressure  pt  are  known.  Then  eq  27  can  be  solved  for  dh/dW,  the  slope  of  the  condition  or 
operating  line  on  the  psychrometric  chart.  A  line  segment  having  this  slope  is  then  drawn  on  the 
chart  through  the  point  representing  the  entering  state  of  the  air.  A  short  distance  along  this  line  a 
new  state  "1  ”  is  arbitrarily  located.  Then,  rewriting  eq  20  as 

=--!0±.(Ah-hf&W)  (28) 

mwCp 
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the  water  temperature  7W  ,  corresponding  to  state  "I  "can  he  calculated.  Nest  eq  27  is  icsolved  lor 
(dhldW  )\  and  the  pioceduic  repeated  until  the  linal  aii  tempeiature  is  determined.  Perlornting  an 
energy  balance  lor  the  entire  ail  column  pi  ovules  a  chec  k  on  the  linal  air  tempeiature  deteimined 
Ironi  the  condition  line  on  the  psyehi omen ic  chart,  or 


ma^ai  =  *  lww-wa(Mu-|V,)l/,l,»- 

Rearranging  yields 

m , 


(29) 


(30) 


It  the  leaving  air  temperature  has  been  determined  and  the  average  mass  transfer  coelficicnt  h^/  1V  is 
known,  then  the  volume  ol  air  mav  be  obtained  from 


m ,  f 

l/=  / 

Mv  l 


dW 
Wr  It' 


(31) 


The  above  integral  can  be  solved  numerically  or  graphically.  (One  convenient  graphical  method  is 
the  use  of  the  Stevens  diagram.) 

For  this  study,  a  computer  program  was  written  (App.  A)  to  calculate  leaving  air  temperature 
and  humidity  ratio,  leaving  water  temperature,  and  percent  moisture  added  based  on  the  loregoing 
analysis  and  the  calculation  procedure  presented  below.  The  input  variables  to  the  computer  pro¬ 
gram  arc  the  entering  air  temperature,  humidity  ratio  and  velocity,  the  entering  water  temperature, 
drop  size  and  velocity,  flow  ratio  and  air  column  height.  The  program  begins  by  assuming  a  leaving 
water  temperature  at  the  bottom  of  the  column.  Then  based  on  the  given  entering  air  conditions 
and  this  assumed  leaving  water  temperature,  the  saturation  humidity,  thermal  diltusivity,  viscous 
diffusivity,  mass  diffusivity,  and  thermal  conductivity  arc  calculated  at  these  conditions  (-50°  C  to 
+20° C)  according  to  the  following  formulas: 


*s  .  O.OO3894e0  09444'7w 

(kgw/kga) 

(32) 

a  =  1/(57736  -  585.7877,) 

(m2/s) 

(33) 

k  =  1/(80711.7  -  766.157,) 

(m3/s) 

(34) 

D  =  2.227x  10-f‘+273  r 
\  273  ' 

(mJ/s) 

(35) 

k  =  0.024577  t  9.027  x  10_s  7,  (W/mC°)  (36) 

These  equations  were  developed  using  a  linear  least-squares  tit  of  property  data  given  in  Karlckar 
(1977)  for  o,  v,  and  k,  and  IVS  values  were  obtained  from  ASHRAE  (1977).  The  mass  diffusivity 
relationship  is  given  in  Eckert  (1959). 

In  order  to  calculate  the  heat  and  mass  transfer  coefficients  for  the  water  drops  required  for  eq 
22  and  27,  the  Prandtl,  Pr,  Schmidt,  Sc,  and  Reynolds,  Re,  numbers  arc  calculated  based  on  the  pre¬ 
viously  defined  property  data.  The  Nusselt  number,  Nu,  for  heat  transfer  and  Sherwood  number, 
Sh,  for  mass  transfer  cannot  be  calculated  according  to  Ranz  (1952),  eq  1 1  and  12.  At  this  point 
A  IF1  across  the  differential  element  can  be  calculated  according  to  eq  22.  Then  using  this  value,  Ah 
for  the  air  stream  can  be  determined  using  eq  27.  Finally,  A  Tw  can  be  found  using  eq  28.  The 
humidity  ratio,  air  stream  enthalpy,  and  water  temperature  arc  then  incremented  to  determine  the 
conditions  on  the  other  side  of  the  differential  element.  These  conditions  become  the  initial  condi¬ 
tions  for  the  second  element  as  the  process  is  continued  step  by  step  downward  until  the  bottom  of 


T«mptrotur«  (*C) 


/  igure  S.  Cooling  curves  tor  a  system  of  water  drops  falling  ftuni 
vanot/s  heights  and  at  various  flow  ratios  (I'.K.)  in  a  -1 8° C environ¬ 
ment. 


2  3  4  5 

flow  Ratio  (kga/kga) 

I  igure  6.  final  temperature  of  water  drops  falling  20 
in  a  -1 8°C  environment  as  a  function  of  flow  ratio. 


ilu-  column  is  reached.  The  accuracy  of  Ihc  conditions  calculated  by  this  procedure  at  the  bottom 
"I  the  column  can  then  be  checked  by  using  eq  30.  If  greater  accuracy  is  desired,  smaller  step  sizes 
(volume  increments)  can  be  chosen  and  the  entire  procedure  repeated. 

Cooling  curves  lor  a  system  ol  water  drops  fallingfrom  various  heights  in  a  -1 8°C  environment 
are  shown  tn  I  igures  S  and  6.  In  figure  3  the  effects  of  falling  distance  and  drop  size  on  the  final 
temperature  ol  the  drops  are  shown  for  a  flow  ratio  of  0.1  (kgw/kga).  As  anticipated,  drops  falling 
a  greater  distance  cool  more.  Also,  as  the  size  of  the  drops  decreases,  the  thci  mal  equilibrium  condition 
between  drops  and  the  air  is  approached  such  that  drops  tailing  20  m  or  60  m  have  the  same  final 
temperature.  I  his  is  due  to  the  tact  that  equilibrium  is  reached  before  the  drops  travel  20  m  so  that 
the  air  temperature  and  drop  temperature  arc  equal  and  no  further  cooling  occurs.  Also,  since  the 
tlrops  are  smaller,  their  time  ol  fall  is  greater,  providing  more  time  for  cooling  to  take  place.  Cool¬ 
ing  curves  lor  How  latios  ol  0.6  and  1.1  are  also  shown  in  this  figure.  As  expected,  when  the  flow 
ratio  increases,  the  cooling  of  the  drops  decreases  since  more  water  is  contained  within  the  air  column. 

I  his  tact  is  ilcarlv  shown  in  f  igure  6  in  which  final  drop  temperature  as  a  function  of  flow  ratio  is 
presented. 
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APPENDIX  A.  FORTRAN  IV  PROGRAM  TO  CALCULATE  FINAL  DROP  TEMPERATURE  AIR 
TEMPERATURE,  AND  HUMIDITY. 

Propam  is  written  in  British  engineering  units. 


Id  REAL  LE 
20  DATA  I  XT / " Y  " / 

Til  104  CONTINUE 

4U  PRINT,"  HEAT  AND  MASH  TRANSFER  FROM  DROPS" 

50  PRINT,  "  " 

50  PRINT,  "AIR  COLUMN  HEIGHT  (FT)" 

7  i?l  READ,  SH 

RiT  PRINT,  "INLET  AIR  TEMPERATURE  (F)" 

90  READ,  TGI 

iau  PRINT,  "INLET  AIR  VELOCITY  (FPS)" 

110  READ,  VG 

12iT  PRINT,  "INLET  HUMIDITY  RATTO  (LBW/LRA)  " 

130  READ,  WI 

140  PRINT,  "DROP  DIAMETER  (IN)" 

150  READ,  D 
160  D  =  D/12. 

170  PRINT,  "SPRAY  VELOCITY  (FPS)" 

180  READ,  VW 

190  PRINT,  "INLET  WATER  TEMPERATURE  (F)" 

200  READ,  TWT 

210  PRINT,  "INITIAL  FLOW  RATIO  {LBW/LBA ) " 

220  READ,  FI 

230  PRINT,  "FINAL  FLOW  RATIO  (LBW/LBA)" 

240  READ,  FF 
250  FF  *  FF+.01 
250  PRINT,  "" 

270  FR  =  FI 
280  70  CONTINUE 

290  20  X  =  SH/S0 

300  TGI  =  TGI 
310  WG1  =  WI 
320  TW1  =  TWI 

330  ND  =  FR*5/(3. 14lS*D**3)/62.  4 

340  AV  =  ND*3.  1416*0*0*39. 715/(TGT+4S0) 

350  ID  =  0 

360  NEND  =  S0 

370  DO  80  7=1,  NEND 

380  CW  =  1.0 

390  TA  =  (TGl+TWl)/2. 

400  Dt  =  0*12. 

420&  0^**4  2886*01-2719.  527*01*01+9164. 955*01**3-12738. 54* 

430  VR  =  VW+VG* (TG 1+4  6  0)/ (TG 1+4  60 ) 

440  VK  =  . 4542439/3600*EXP ( .0039815*TA) 

4S0  TO  =  (.63571 7R+.0025373*TA)/3600 
460  PR  =  VK/TD 
470  RE  =  VR*D/VK 

480  TK  =  ( .013115  +  2.63929E-5*TA)/3600 
490  HC  =  TK*(2.0+.6*PR**.333*SORT  (RE))/0 
S00  OEN  =  39 . 71S/ (TG 1+460 ) 

510  IF(IO.EO.l)  GO  TO  102 

520  101  OC  =  (2. 397E-4)*( (TA+460)/4B0)**i.ri 

530  GO  TO  103 

540  102  OC  »  (3.86779F.-9)*(  (TA+460)*B./q.)**i.941s 
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5S0  103  CONTINUE 

860  sc  *  vk/oc 

S70  CP  *  .24+. 45*WG1 

Sft«  HO  *  DC* (2. 0+. 6*SC** . 333*SORT (RE) ) /O 
890  LE  *  HC/ (HD*DEN*CP) 

800  WS  =  .  000726S*EXP  ( .0824685*TW1) 

610  HWG  =  107S+.  4S*TG  1 
620  HG1  =  ( .24+. 45*WG1) *TG1+1075*WG1 
^30  HGSW  *  ( .240  +  . 45*WS) *TW1+107S*WS 
<540  OWG  =  DEN*HD*AV*X*  (WS-WG1) 

660  HFW2  «  CW*  (TW2-32  ) 

660  OHG  =  DWG*  (LE*  (HG 1-HGSW)  / ( WG1-WS )  +HWG-1075*LE ) 
6™  OTW  =  +DWG/(FR*CW) *HFW2-DHG/(FR*CW) 

680  TW2  =  TW 1  +  OTW 

690  WG2  =  WG1  +  OWG 

700  HG  2  =  HG1  +  OHG 

710  TG  2  =  (HG 2-1 07 S*WG2) / ( . 24  +  . 45*WG2) 

720  TGI  =  TG  2 

730  WG1  =  WG2 

740  TW1  =  TW2 

760  80  CONTINUE 

760  o  *  0*12. 

770  PRINT,  "OROP  DIAMETER  *",  D 

780  PRINT,  "FLOW  RATIO  »",  FR 

790  PRINT,  "FINAL  WATER  TEMPERATURE  (F)",  TW2 

800  PRINT,  "LEAVING  AIR  TEMPERATURE  (F)",  TG2 

810  PRINT,  "FINAL  HUMIOITY  RATIO  (LBW/LBA) ",  WG2 

820  PRD  =«  (  (WG2-WH/WI)*100. 

830  PRINT,  "PERCENT  WATER  AOOED" ,  PRO 
840  PRINT,"  " 

8S0  PR  =  FR+.S 
860  0  *  0/12. 

870  IF  (FR.LT.FF)  GO  TO  70 

880  PRINT, "NEW  DATA  "*  YES-Y  OR  NO-N" 

890  READ (8 , 105 )  TANG 
900  108  FORMAT ( A1 } 

910  IF (IANS. EO. IXT)  GO  TO  104 
920  STOP 
930  END 
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A  facsimile  catalog  card  in  Library  of  Congress  MARC 
format  is  reproduced  below. 


Zarling,  John  P. 

Heat  and  mass  transfer  from  freely  falling  drops 
at  low  temperatures  /  by  John  P.  Zarling.  Hanover, 
N.H.:  U.S.  Cold  Regions  Research  and  Engineering 
Laboratory;  Springfield,  Va. :  available  from  Nation¬ 
al  Technical  Information  Service,  1980. 

v,  20  p.,  illus.;  28  cm.  (  CRREL  Report  80-18.  ) 
Prepared  for  Directorate  of  Military  Programs, 
Office,  Chief  of  Engineers  by  Corps  of  Engineers, 
U.S.  Army  Cold  Regions  Research  and  Engineering 
Laboratory  under  DA  Project  4A762730AT42 . 
Bibliography:  p.  10. 


Zarling,  John  P. 

Heat  and  mass  transfer  from  freely  falling... 

1980  (Card  2) 

1.  Drops.  2.  Freezing.  3.  Heat  transfer.  4.  Ice. 
5.  Mass  transfer.  I.  United  States.  Army.  Corps  of 
Engineers.  II.  Cold  Regions  Research  and  Engineering 
Laboratory,  Hanover,  N.H.  III.  Series:  CRREL  Report 
80-18. 
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